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\Q ■ We report here results on charm and strange hadron spectra. Adopting a mixed action approach, 

we use overlap fermions for valence quarks, on a background of 2+1+1 flavours HISQ gauge 
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1. Introduction 

During the last few years MILC lattice collaboration has generated [JTJ] a large set of con- 
figurations with the one-loop, tadpole improved Symanzik gauge action and the highly improved 
staggered quark (HISQ) fermion action [Q]. Four flavours of dynamical sea quarks: degenerate up 
and down, strange, and charm, were included in these configurations. Taste violations in the HISQ 
action were found to be small [Qj. On the other hand, overlap fermion action []3|] has exact chiral 
symmetry [||, |J] on the lattice and is automatically &{a) improved. It has also many desirable 
features, such as the adaptation of multi mass algorithms []|]. However, using overlap action for 
the dynamical quarks is still prohibitively costly, except with fixed topology ||J|. In this work we 
have adopted a mixed action formalism where we used overlap valence quarks on the HISQ gauge 
configurations generated by MILC collaboration [p. By adopting such a mixed action approach, 
one can get advantage of the chiral symmetry and low quark mass limit of overlap fermions, and the 
advantage of having a large set of configurations with small discretization errors as well as small 
taste breaking effects. One also gets the advantage of simulating both light, strange as well as heavy 
fermions on the same lattice formalism with chiral fermions having no (a) errors. Of course we 
will have to encounter the usual issues related to mixed action and partial quenching Jl4[]. Similar 
approach has been taken by the %QCD collaboration, who used overlap valence quarks on 2+1 
flavours dynamical domain wall gauge configurations [0]. 

In this work, by using above mentioned mixed action approach, we report our preliminary 
results on charm and strange hadron spectra as well as leptonic decay constants for D s and D* 
mesons. 

2. Numerical details 

We used two sets of dynamical 2+1+1 flavours HISQ lattice ensembles, generated by the 
MILC collaboration : a set of 32 x 96 lattices with lattice spacing a = 0.09 fm, and a set of 
48 3 x 144 lattices with a = 0.06 fm. The details of these configurations are summarized in Ref. |jl|]. 
The results reported here were obtained from 100 configurations on the coarser lattice, and 35 
configurations on the finer lattice. Here we did not address the uncertainties in the determination 
of lattice spacings which we will address in subsequent works. 

For valence quarks we used overlap action [Q]. For the numerical implementation of massive 
overlap fermions the methods used by the %QCD collaboration [JSj] were followed. To evaluate the 
sign function, the Zolotarev approximation is used, and the low Wilson eigenmodes are projected 
out by Arnoldi method. We used 350 and 160 Wilson eigenmodes for the coarser and the finer 
lattices respectively. The constraint here came from the size of the memory of the machine. Using 
multimass CG algorithm we were able to input a range of quark masses covering from light, strange 
to heavy quarks. It is worthwhile to mention that one should not use the same precision for the 
evaluation of light and heavy quark propagators. While light quark and strange quark propagators 
needs relatively lower precision, for heavy quark propagators which decay rapidly, having higher 
precision in the inversion is absolutely necessary. Otherwise, one will get flat hadron correlators at 
the middle of the lattice. However, for using different precisions in the multimass we found that 
it is not necessary to reload the eigenmodes repeatedly, which increases the total inversion time. 
We used periodic boundary condition in the spatial and antiperiodic in the temporal directions. 
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Gauge configurations were first fixed to coulomb gauge and then HYP smeared (one HYP). Using 
both point and wall sources we calculated various point-point, wall-point as well as wall-wall 
correlators. 

The strange mass was tuned by setting the ss pseudoscalar mass to 685 MeV [^]. With the 
strange mass tuned this way, the mass of the vector <p(ss) meson was found to be 1050(20) MeV 
and 1030(10) MeV on the coarser and finer lattices respectively, and the mass of the £l(sss) baryon 
was found to be 1620(40) and 1630(30) MeV respectively. The charm mass is tuned by setting 
the spin-averaged IS state mass, (m^ c + 3mj/ w )/A, to its physical value. The charm quark mass 
{ma c harm) hi these lattices are obtained to be 0.55 and 0.336, respectively. Interestingly, the MS 
mass of charm, 1.275 GeV, corresponds to ma = 0.582 and 0.388, respectively, indicating that the 
mass renormalization constant Z^ s is close to 1 for this action. 

Since ma is not very small, we need to be careful about discretization errors. Overlap action 
does not have 0{ma) errors. In order to estimate the size of discretization errors coming from 
higher orders of ma, we look at the energy-momentum dispersion relation of the IS charmonia. We 
calculated the mesons at various external momenta p 2 = (2n/L) 2 n 2 , with n < 2. We used both point 
source and a momentum-induced wall source for the finite momenta mesons. For the wall source at 
finite momenta, we put a suitable phase factor in the wall so as to project to the suitable momentum. 
While this requires a separate inversion for each momenta, that is more than compensated by the 
improvement in the signal, as shown in Fig. 1. In both figures we show both point source and wall 
source results for the coarser lattice. Note that the wall source results are from 8 configurations, 
but the statistical error is already comparable to that from 100 configurations with the point source. 

In the left panel of Fig. 1, we show E(p) for various momenta for the pseudoscalar meson 
on our coarser lattices. Also shown are the continuum dispersion relation, E 2 = m 2 + p 2 , and 
the lattice dispersion relation for the standard scalar action with 0{ma) 2 error, sinh 2 (E(p)a/2) = 
(sinh 2 (ma , /2)) 2 + sin 2 (pa/2). While substantial deviation from the continuum dispersion relation 
is seen, the lattice scalar dispersion relation seems to explain the data quite well. For a more 
quantitative analysis, we follow the standard practice of introducing an effective "speed-of-light" c 
through E 2 (p) = m 2 + p 2 c 2 . The value of c obtained using this relation and E(p) at various p are 
shown in the right panel of Fig. 1. As the figure shows, we get c 2 ~ 0.75 on our coarser lattices, and 
c 2 ~ 0.89 on the finer lattices. Using quenched overlap fermions, on an even coarser lattice, Ref. 
[10] found c 2 ~ 1 for the charm. However, on our dynamical lattices, we obtained c 2 values similar 
to that obtained in the literature with clover action, indicating a similar size of cutoff in the overlap 
action as that in the clover. The c 2 values are well-approximated by the estimate ma/ sinh(rafl r ), 
suggested by the lattice scalar action. 



3. Results 

We utilized both point and wall sources to generate a large set of quark propagators. Us- 
ing those we calculated various hadron correlators and extracted light, strange and charm hadron 
masses. The pseudoscalar meson masses ranged between 240 — 3582 MeV for the coarser lattices 
and between 420 — 5275 MeV for the finer lattices. In Fig. 2 we plot the pseudoscalar meson 
masses over a range of quark masses. 
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Figure 1: (Left) Energy-momentum dispersion relation for the tj c on the lattices with a = 0.09 fm. Also 
shown are the continuum relativistic dispersion relation (solid line), and the lattice dispersion relation 
(dashed line) for the standard scalar action. (Right) c 2 = (E 2 (p) —E 2 (p = 0))/p 2 calculated for r\ c at 
various values of p, on lattices with a = 0.09 and 0.06 fm. Solid horizontal lines represent the estimate 
ma I smh(ma), suggested by the lattice scalar action. 



Using these correlators, we report here preliminary results for spectra of charmonia, mesons 
of D s family, and charmed baryons. We also report preliminary studies of leptonic decay constants 
of D s and D* mesons. 
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Figure 2: Pseudoscalar meson masses over a range of light to heavy quark masses. Left figure is for m 2 a as 
function of quark mass ma, and the right figure is for m K a vs ma in physical unit. 



3.1 Hyperfine splitting in IS charmonia 

The hyperfine splitting in IS charmonia is one of the most well studied physical quantities 
in lattice charmonium calculations over the year, and until very recently [ JITl ] lattice results were 
found to be smaller than the experimental value (~ 116 MeV). This underestimation was mainly 
due to the discretization error associated with the charm quark action and the quenched approxi- 
mation. In our study we also calculated this splitting. In Fig. 3, we showed the effective splittings 
between vector and pseudoscalar correlators (jackknifed) at the tuned charm mass for wall-point 
correlators. Horizontal lines shown are the fit results obtained by fitting only these correlators. Our 
final estimated results, by simultaneously fitting the wall-point and point-point correlators, for this 
hyperfine splitting are 114+| MeV and 109^3 MeV corresponding to lattices with spacings a = 
0.09 and 0.06 fm respectively. 
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Figure 3: Effective hyperfine splitting in IS charmonia for wall-point correlators for lattices with spacing 
0.09 (left) and 0.06 fm (right). Horizontal lines show the fit results obtained by fitting only these correlators. 

3.2 Charm and strange mesons 

Along with pseudoscalar and vector mesons we also calculated the ground state spectra of 
tensor, axial and scalar charmonia. Besides charmonia we also calculated the strange-charm D s 
mesons. In Fig. 4(a) we showed splittings between these various mesons at two lattice spacings. It 
is to be noted that the splitting of D s — T] c /2, which has smaller lattice spacing error [12], obtained 
is consistent with its experimental value. 



500 

a 

•g 400 



300 









h c -J\ c 
3CcJ-1c • 


t 


V 


Xc0-% • 




Dj-d, • 


\ 


\ ' 


D s - D , 

m - n c • 


• 


• 



> 



3600 



0.02 0.04 

a ( fernii) 
(a) 



d 3400 

s 

"§ 3200 
G 3000 

s 

2800 
2600 







Q. cc baryons (ccs) 


* 






This work 








Briceno et al 








Mathur et al 








Expt 








Ci c baryons (ccs) 














* 


f 






1/2 + 


3/2 + 


1/2" 



(b) 



Figure 4: (a) Meson mass splitting for charmonia and strange-charm mesons at two lattice spacings. Ex- 
perimental values are shown in the left side, (b) Q. c and Q. cc baryon masses. Also shown are other lattice 
determinations, and the experimental values, where available. 

3.3 Charmed baryons 

Recently, there have been exciting developments in heavy-baryon physics, both theoretically 
and experimentally. Results from LHC and future charm-bottom factories are expected to add to 
the excitement in this field in the near future. On our lattices we extracted ground state spectra of 
charmed baryons with one or more charm quark content, for example, baryons with quark content 
csu,cuu,css and ccs. In fig. 4(b) we showed results for Cl c (css) and £l cc (ccs) baryons. It is 
to be noted that for these baryons we extracted masses for both spin 1/2 and spin 3/2 with both 
parities, some of which are yet to be measured experimentally. Our results are consistent with 
other lattice results [O], and the experimental values, where available. Results for other charm 
baryons at various pion masses were shown in Fig. 5. The chiral extrapolation will involve effects 
from partial quenching as well as mixed action formalism. Mixed action partially quenched chiral 



perturbation theory (MAPQ^PT) has been developed for various mixed action formalisms [14] 
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which requires calculation of A„„ T [15], the low energy constant representing Gia 1 ) discretization 



dependence. In future, utilizing MAPQ^PT we will extrapolate these masses to the physical limit. 
We also calculated the experimentally unknown triply charmed baryon £l ccc , and a better way to 
quote its mass is through the splitting of ma ccc — |m//>p, which we extracted to be 1 IO^q MeV and 
120(10) MeV, for 0.09 fm and 0.06 fm lattices respectively. 
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Figure 5: Ground state spectra for charmed baryons at various pion masses for 32 3 x 96, a = 0.09 fm lattice. 
3.4 Decay constants 

The decay constants of heavy-light mesons are experimentally very important: they are es- 
sential ingredients in extracting CKM matrix elements from decays of heavy-light mesons. For 
the suitable continuum currents {A^, V^} = {Za (sYuTfc c )hti Zy (sYu c )ht}i where Za,Zv are the 
renormalization constants to match the lattice operators to continuum currents, the leptonic decay 
constants for the pseudoscalar D s and the vector D* are defined as 

(0\A^\D s (p)) = if D , Pll , <0|V M |D;(p,A)) = ifn, in,,, r^. (3.1) 

where is the polarization vector for the vector meson. From the amplitude of the (V(t) V(0)) 
coiTelator, fjy> can therefore be extracted using standard techniques. Similarly, fj) s can be extracted 
from (A$(t)P{p)) and (P(t)P(O)) con-elators, where P = Zpsy=,c is the pseudoscalar current. 

Employing point-point as well as wall-point and wall-wall correlators we extracted the suitable 
amplitudes for the above matrix elements. The relevant renormalization constants, Zy and Za have 
not been calculated yet. On the other hand, since the overlap action possesses exact chiral symmetry 
on the lattice, the ratio Zy /Za is expected to be close to 1 (it is exactly 1 for massless quarks). This 
allows us to get an estimate of the size of frr through the calculation of 

pff f \ - Za f°* ~ f°t 

K{Jd' s Jd s ) - ^-- t 7— ■ 

Zy Jd s ]D S 

since fo s has been calculated very precisely in the literature. Our preliminary results for R(fo* , /d s ) = 
1.15(10) on the lattices with a = 0.09/;n. The mixed action effects will be smaller for heavy-light 
mesons, and in the ratio its effects will be minimal. 

4. Conclusions 



In this work we reported preliminary results on the ground state hadron masses along with 
charm-strange meson decay constants by using a mixed action approach, comprising overlap va- 
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lence quarks generated on the background of dynamical 2+1+1 flavours HISQ configurations. The 
results, in particular the hyperfine splitting of IS charmonia, are encouraging and suggest that the 
overlap valence on 2+1+1 flavor HISQ configurations is a promising approach to do lattice QCD 
simulation with light, strange and charm quark together in the same lattice formulation. However, 
at these lattice spacings, the discretization errors of the overlap action for the charm, as evident 
from the dispersion relation, are not negligible, and are similar in size to those of the clover action. 

This is a continuing study and we expect to be able to do suitable chiral and continuum extrap- 
olations, to make experimentally relevant predictions for various charmed baryons. The splitting 
(m^ — |m//>p), between ■//*}/ and the unknown triply charmed baryon £l ccc was found to be 
110+^ MeV and 120(10) MeV, on our lattices with a = 0.09 and 0.06 fm respectively. We are also 
studying heavy-light decay constants. Necessary renormalization constants will be calculated in 
future. Our preliminary results for fy* /fo s (for which the renormalization constants approximately 
cancel) is 1.15(10), from our coarser lattice. 

5. Acknowledgement 

The computations were carried on the Blue Gene P of Indian Lattice Gauge Theory Initiative, 
Tata Institute of Fundamental Research (TIFR), Mumbai, and on the Konark cluster, NISER. We 
would like to thank A. Salve and K. Ghadiali for technical support. We are grateful to the MILC 
collaboration and in particular to S. Gottlieb, for providing us with the HISQ lattices. We also 
thank A. Li and A. Alexandru for help with numerical issues of overlap fermions. 

References 

[1] A. Bazavov et al. Phys. Rev. D 82, 074501 (2010); A. Bazavov et al, PoS(Lattice 2010)320 (2010); 
A. Bazavov et al. (MILC Collaboration), PoS(Lattice2012)158. 

[2] E. Follana et al, Nucl. Phys. B (Proc. Suppl.) 129 and 130, 447 (2004); E. Follana et al, Phys. Rev. 
D75, 054502 (2007). 

[3] H. Neuberger, Phys. Lett. B417 (1998) 141; ibid. B427 (19998) 353. 

[4] M. Luscher, Phys. Lett. B428 (1998) 342. 

[5] R. Edwards et al, Phys.Rev. D59 (1999) 094510. 

[6] H. Fukaya et al, [JLQCD Collaboration], Phys. Rev. Lett. 98, 172001 (2007); Phys. Rev. D 77, 
074503 (2008). 

[7] A. Li et al Phys. Rev. D82 (2010) 1 14501; N. Mathur et al, PoS LATTICE2010 (2010) 1 14. 

[8] Y. Chen et al, Phys. Rev. D70 (2004) 034502, S.J. Dong et. al, Phys. Rev. Lett. 85 (2000) 5051-5054, 
ibid. Phys.Rev. D65 (2002) 054507. 

[9] C.T.H. Davies et al, Phys.Rev. D81 (2010) 034506. 

[10] S. Tamhankar et al, Phys. Lett. B638 (2006). 

[11] T. Burch et al, Phys.Rev. D81 (2010) 034508; G.C. Donald et al, Phys.Rev. D86 (2012) 094501. 
[12] C. Davies, et al, Phys.Rev. D82, 1 14504 (2010). 

[13] N. Mathur et al, Phys.Rev. D66 (2002) 014502; R. Briceno et al, Phys.Rev. D86 (2012) 094504. 

[14] O. Bar et al, Phys. Rev. D 67, 1 14505 (2003), ibid. 70, 034508 (2004); J.W. Chen et al, Phys. Rev. 
D75. 054501 (2007); K. Orginos et al, Phys. Rev. D77, 094505 (2008) and references therein. 

[15] M. Lujan et al, Phys.Rev. D86 (2012) 014501. 



7 



